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Abstract. Zeros of the nth moment of the partition function [Z"] are investigated in 
a vanishing temperature limit (3 oo, n keeping y = [3n ~ 0{1). In this limit, the 
moment parameterized by y characterizes the distribution of the ground-state energy. 
We numerically investigate the zeros for ±J Ising spin glass models with tree and 
other several systems, which can be carried out with a feasible computational cost by 
a symbolic operation based on the Bethe-Peierls method. For several tree systems we 
find that the zeros tend to approach the real axis of y in the thermodynamic limit 
implying that the moment cannot be described by a single analytic function of y as 
the system size tends to infinity, which may be associated with breaking of the replica 
symmetry. However, examination of the analytical properties of the moment function 
and assessment of the spin-glass susceptibility indicate that the breaking of analyticity 
is relevant to neither one-step or full replica symmetry breaking. 
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1. Introduction 

Spin glasses are a typical example of disordered systems and have been investigated for 
a long time The first comprehensive understanding of spin glasses was obtained 
by investigating the so-called SK model introduced by Sherrington and Kirkpatrick 
[2], which describes fully connected Ising spin glasses. In analyzing this model, they 
employed the replica method under the replica symmetric (RS) ansatz. However, the 
SK solution contains an inconsistency in that the entropy at low temperatures becomes 
negative. This problem has led to much controversy regarding the validity of the replica 
method. In 1980, Parisi developed the replica symmetry breaking (RSB) scheme [Sj H] 
and showed that a sufficient solution can be obtained within the framework of the replica 
method. 

Although Parisi's RSB scheme is consistent at low temperatures, a mathematical 
justification of the replica method and a proof of the Parisi scheme were lacking until a 
recent study showed that the Parisi's solution is exact for the SK model [5|. However, 
this does not resolve all of the questions regarding the replica method. There are still 
many unsolved issues, e.g. ultrametricity and the origin of the RSB. These issues have 
attracted renewed interest as applications of the the replica method have increased 
rapidly O [8], |9], and a deeper understanding of this method is greatly desired. 

The RSB is considered to relate to the analyticity of a generating function g{n) 
defined as follows: 

g{n) = lim ^jv(^), (1) 

N^oo 

g^(n)^j^log[Z% (2) 

where n is referred to as the replica number and the brackets [■ ■ ■] denote the average 
over the quenched randomness. The functions gN{n) and g{n) are defined for Vn e M 
(or G C) and the free energy is derived from gN{n) as 

f = — lim [logZl = — lim lim — QNin). (3) 

The name 'replica method' is often used to indicate the second identity, though this 
method should be considered as a systematic procedure to evaluate eqs. ^ and ([2]). 
In general, the calculation of [Z^] is difficult for real n G M (or complex C). To 
overcome this difficulty, the replica method first computes [Z"] for natural numbers 
n = 1,2, ■■■ G N, then extends the obtained expressions of [Z"] to n G M by their 
analytical continuation. However, this technique causes the following two problems. The 
first concerns the uniqueness of the analytical continuation from natural to real numbers. 
Even if all the moments of [Z"] are given for n G N, in general it is impossible to uniquely 
continue the analytical expressions for n G N to n G M (or C). Carlson's theorem 
guarantees that the analytical continuation from n G N to n G C is uniquely determined 
if [Z""]^^^ < ©(e'^'"') holds as Re(n) tends to infinity [lOj. Unfortunately, the moments 
of the SK model grow as 

[^n]i/iv < o(e^l"l'), where C is a constant, and therefore 
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this sufficient condition is not satisfied, van Hemmen and Palmer conjectured that the 
failure of the RS solution of the SK model might be related to this issue, though further 
exploration in this direction is technically difficult pJ!]. The second issue concerns the 
possible breaking of the analyticity of g{n). In general, even if gN{n) is guaranteed to be 
analytic with respect to n for finite N, the analyticity of g{n) = limTv^oo fi'Af('T-) can be 
broken. Since it is unfeasible to exactly compute [Z"] except for a few solvable models, 
in most cases, only the asymptotic behavior is investigated by using certain techniques 
such as the saddle-point method in the limit — * oo. This implies that, in such cases, 
the expression analytically continued from n G N to n G M in the limit ^ oo will 
lead to an incorrect solution for n ^ if the breaking of analyticity occurs in the 
region < n < 1. Recently, it has been shown that analyticity breaking does occur and 
is relevant to one-step RSB (IRSB) for a variation of discrete random energy models 
[T2l [T3l [131 US] , for which the uniqueness of the analytical continuation is guaranteed by 
Carlson's theorem and for which [Z^] or equivalently gN{n) can be assessed in a feasible 
manner without using the replica method for finite and n G C. This is a strong 
motivation to investigate the analyticity of [Z"] for various systems to explore possible 
links to different types of RSB. 

Under this motivation, we investigate possible scenarios of analyticity breaking of 
g{n) = limiy^oo gN{n). For this purpose, we observe the zeros of [Z"-], which will be 
referred to as "replica zeros" (RZs), on the complex plane n G C for finite and 
examine how some sequences of zeros approach the real axis as A^ tends to infinity. 

For the discrete random energy model mentioned above, this strategy successfully 
characterizes an RSB accompanied by a singularity of a large deviation rate function 
with respect to A^^^ logZ [16j. As other tractable example systems, we investigate ±J 
models with a symmetric distribution on two types of lattices, ladder systems and Cayley 
trees (CTs) with random fields on the boundary. There are two reasons for using these 
models: Firstly, these models can be investigated in a feasible computational time by the 
Bethe-Peierls (BP) approach [17]. Especially, at zero temperature this approach gives a 
simple iterative formula to yield the partition function. Employing the replica method 
and the BP formula, we can perform symbolic calculations of the replicated partition 
function [Z^], which enables us to directly solve the equation of the RZs [Z"] = 0. The 
second reason is the existence of the spin-glass phase. It is known that the spin-glass 
phase is present for CTs |T8l [191 ISQl [21] and is absent for ladder systems. Therefore, 
we can compare the behavior of RZs, which are considered to be dependent on the 
spin-glass ordering. 

This paper consists of five sections. In the next section, we give an explanation 
of our formalism. Simple recursive equations to calculate [Z"] are derived in a zero- 
temperature limit by combining the BP approach and the replica method. The 
relationships of CTs to Bethe lattices (BLs) and regular random graphs (RRGs) are 
also argued. Assessing the contribution from the boundary indicates that IRSB does 
not occur in CTs and BLs while it does for RRGs in the thermodynamic limit when the 
boundary contribution is correctly taken into account, which is the case in the evaluation 
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of RZs. This implies that the possible RZs of a CT are irrelevant to IRSB. In sec. 3, 
we present plots of RZs and investigate their behavior. Their physical significance is 
also discussed. In sec. 4, a possible link to another type of RSB, the full RSB (FRSB), 
is examined. Numerical assessment of the de Almeida-Thouless (AT) condition based 
on the divergence of spin-glass susceptibility, however, indicates that RZs do not reflect 
FRSB, either. Therefore, we conclude that the analyticity breaking that occurs in CTs 
is irrelevant to RSB. The final section is devoted to a summary. 



2. Formulation 



In this section, the main ideas of the paper are presented. It is shown that the RZ 
equation [Z^] = is simplified at zero temperature. An algorithm to evaluate the 
generalized moment [Z^] for n G C is developed by introducing the replica method to 
the BP approach. 



2.1. Equation of the replica zeros at zero temperature 
Solving 

in = (4) 

with respect to n is our main objective. Unfortunately, this is, in general, a hard 
task even by numerical methods because eq. (jl]) is transcendental and becomes highly 
complicated as the system size grows. In the T — limit, however, the main 
contributions to the partition function come from the ground state and eq. (jlj) becomes 

[Z"] ^ [rf^e"^"^-] = 0, (5) 

where Eg is the energy of the ground state and dg is the degeneracy. If n is finite 
when /3 — s> oo, the term e"^"^*^ diverges or vanishes and there is no meaningful result. 
Therefore, we suppose that non-trivial solutions exist only in the limit n — 0, /? — oo, 
and y = (3n ~ 0(1). This assumption is consistent with the fact that the solution of 
the SK model is well-defined in this limit [22]. Under this condition, eq. ([5]) becomes 

[e-^^^] = 0. (6) 

In the following, we focus on the ±J model whose Hamiltonian is given by 

if = - ^ JijSiSj, (7) 

and the distribution of interactions is 

P{J^J) = l6iJ,,-l) + hiJ,, + l), (8) 

assuming that the total number, Nb, of interacting spin pairs {i,j) is proportional to A^, 
which is the case for ladder systems and CTs. This limitation restricts the energy of any 
state to an integer value. As a result, eq. can always be expressed as a polynomial 
of a; = e^, which significantly reduces the numerical cost for searching for RZs. 
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One issue may be noteworthy here. In the present study, we focus on the hmit 
ri — > 0, /5 — > CX3 keeping i3n ^ y ~ 0{1). In research on zeros of partition functions, on 
the other hand, another limit n keeping (3 finite can be examined as well. In the 
latter case, the zeros with respect to complex f3 are sometimes referred to as "Fisher 
zeros" [23]. Intuitively, Fisher zeros characterize the origin of singularities with respect 
to f3 for typical single sample systems. These can be examined not only for random 
systems [T3l [H] but also for systems of deterministic interactions such as frustrated 
anti- ferromagnetic Ising spin models [21]. As n — limit is taken on ahead for each 
/3, Fisher zeros are irrelevant to the analyticity concerning the replica number n. For 
examination of the analyticity with respect to n, it is necessary to investigate the zeros of 
[Z"] in the complex plane of n. In the situation of vanishing temperatures jS —>■ oo, this 
naturally leads to the current nontrivial limit n — > 0, /? — oo keeping f3n y ~ 0(1). 



2.2. The Bethe-Peierls approach 

2.2.1. General formula On cycle free graphs, it is possible to assess the partition 
function by an iterative method, i.e. the BP approach. We here present a brief review of 



the procedure for CTs. The BP approach in ladder systems is presented in Appendix A 
The basis for our analysis is a formula for evaluating an effective field by a partial 
trace: 

"^ex^ {(^{J'ijSiSj + hjSj)] = Aey.Y>{(^hiSi). (9) 

So 

A simple algebra offers 

^^2cosh/jJ,.cosh/?/^,^ ^^^^ 



where 



cosh (3hj 



(3h j = ianh. ""^(tanh/^Jjj tanh/J/ij). (11) 



The fields hj and hj are sometimes termed the cavity field and cavity bias, respectively. 
For CTs, iterating the above equations from the boundary gives the series of cavity 
fields and biases {hj, hj}. In general, a cavity field becomes a summation of the cavity 
biases from its c — 1 descendants (c is the coordination number): 

c-l 

h. = Y.h,. (12) 

i=i 

Hereafter, we mainly focus on the c = 3 C3iSG, clS shown in fig. [H but the extension 
to general coordination numbers is straightforward. In addition, generalizing to fc-spin 
interacting CTs (fc-CTs) is also straightforward; the only necessity is to replace the 
partial trace ([9]) with that for a fc-spin interaction, as 

exp U [SiJk \{S, + Y. ^^^A I = AeM^h^S,), (13) 

Si,52 I \ 3=1 3=1 J ) 
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2^-^ cosh /5 Jfc Yl)-] cosh ph. 

K = h,, A= ^^r^ (14) 

cosh (3hk 

hk = — tanh ^ j tanh/JJ^ JJ^ tanh/5/ij j . (15) 




Figure 1. Local structure of a CT with coordination number c — 3. 

Let us denote the partition function in the absence of i's ascendants as Zj. 
Equations (P)- (IT^ imply that the partition function is updated as 

= ^ ZjZk exp{-P{AHij + AHik)}pj{Sj, hj)pk{Sk, hf,), (16) 

Si^Sj 

where 

2cosh/3/., ' ^^^^ 
is the one-site marginal in the absence of j's ascendants and Aifjj = —JijSiSj is the 
bond Hamiltonian added by a propagation procedure. 

As a final step, the contribution from the origin of the tree is calculated as 

3 

Z = Z1Z2Z3 JJ(2cosh/?Ji) 

i=l 

/ 1 + tanh/3 Ji tanh/9 J2 tanh jShi tanh l3h2 + (two terms with 1, 2, 3 rotated) 

X I 

V 4 
and the whole partition function Z is derived. Taking the T — > limit yields the 
ground-state energy for a given bond configuration. In this limit, eqs. ( |TT|) and ( |T6l) 
become 

hj sgn{Jijhj) mm{\ Jij\, \ hj\), (19) 

and 

lim --logZi = Ei 

_ p , p 7 7 , / ( sgn{JijJikhjhk) > ) 

I 2 mm( I Jij \,\Jik\,\hj\,\nk\) [ otherwise ) 

We assume sgn(O) = in this paper. 
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2.2.2. Use of the replica method For a given single sample of interactions and boundary 
conditions, a simple application of the BP algorithm enables us to evaluate the partition 
function in a feasible computational time. Unfortunately, this does not fully resolve the 
problem of the computational cost for assessing the moments (jlj) since the cost for 
evaluating an average over all possible samples of interactions and boundary conditions 
grows exponentially with respect to the number of spins. However, this difficulty can be 
overcome by analytically assessing the configurational average for n G N and analytically 
continuing the obtained expressions to n G C in the level of the algorithm, a method 
which may be considered a generalization of the replica method. 

For this purpose, we first evaluate the nth moment of the partition function Z 

n 

E{n) ^ [Z-] = Tr Hiexp^, S^Sf)], (21) 

for n G N, where a is the replica index. Let us denote the effective Hamiltonian as 

i/eff = J2 = 1^-^ log[exp(/?J,, ^ SrS^)], (22) 

where [■ ■ ■] stands for the configurational average with respect to the interactions {Jij}. 
This means that eq. (I2T1) is simply the partition function of an n-replicated system, 
which is defined on a cycle free graph and is free from quenched randomness. Therefore, 
by expressing the BP algorithm in the current case as 




S,(n) = Yl 

S i,S j ,S k 
= Ypi{Si)Ei{n 

where pi is the one-site marginal distribution of site i, eq. (1211) can be assessed in a 
feasible time. The expressions f l23|) and fl24|) define the updating rules of pi and Sj(?T,). 

So far, we have made no assumptions or approximations and therefore eq. fl24l) 
yields exact assessments for n G N, given a boundary condition. To generalize this 
scheme to n G C, we here introduce the RS ansatz, which is the second step of the 
replica method and, in general, is expressed by a restriction of the functional form of 
Pi{Si) as 

ra^ f TT A + tanh(/?/i)5f \ ^ f e'"'^'^^? 
P.iSd = j Mh) n [ ^ )dh = j -mj^^^^^^Jh, (25) 

where 7rj(/i) is a distribution to be updated in the algorithm. The expression of eq. 
(!25l) guarantees that Pi{Si) is invariant under any permutation of the replica indices 
a = 1, 2, . . . , n. Note that this property is automatically satisfied over all of the objective 
lattice if only the distributions on the boundary are expressed in the form of eq. fl25|) . 
Inserting eq. (l25l) into eq. ( 12^ and performing some simple algebraic steps gives 

/f3hi 5" f r r 
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(26) 
(27) 



(2cosh/5/ii)" 

Equations fl26|) and fl27|) provide an expression of the rephca symmetric BP algorithm 

2 cosh phi 



T^i{hi) (X T^j{hj)TTk{hk) 



6{hi - hj - hk) 



2 cosh i3hj2 cosh jShk 

2 cosh (3{hj + hk) 



dhjdhk,{28) 



(29) 



^jSfc(2cosh^)^" / / dhjdhk7rj{hj)7rk{hk) ^ 

2 cosh i3hj2 cosh jShk 

which is apphcable to Vn G C When the algorithm reaches the origin of the CT, the 
moment of eq. (l2T]) is assessed as 



E{n) = [ZT =SiS2S3(2 cosh /?)3"yyy dhidh2dhs7T,{h,)7T2{h2)Mh) 

/ 1 + tanh/3Jitanh/3J2tanh/3/iitanh/3/;,2 + -R\" , . 

X ^ j J , (<^0) 

where R is two terms with the indices 1, 2, 3 rotated. 

2.2.3. Zero-temperature limit Under appropriate boundary conditions, the zero- 
temperature limit (3 —>■ oo, n — > keeping y = (3n finite, which we focus on in the 
present paper, yields further simplified expressions of the BP algorithm. For this, we 
generate replicated spins of each site on the boundary with an identical random external 
field hi = ±1, the sign of which is determined with an equal probability of 1/2. This 
yields the cavity field distribution 



n,{h,) = -{6{h,-l)+6{h, + l)) 



(31) 



and the partition function 



= (2cosh/5)" ^ e^ (32) 

as the boundary condition. The relevance of the boundary condition to the current 
objective systems is discussed later. 

Equation (I3T1) in conjunction with the property \Jij\ = 1 allows TXi{hi) in eq. (l28ll 
to be expressed without loss of generality as 

c-l 

T^i{hi) = Pr,oS{hi) +Y,Pr,f {^h, -f)+ 5{h, + /)) , (33) 

/=1 

where Pj = {pi-Q,Pi-i, ■ ■ ■ ,Pi-c-i) represents a probability vector satisfying pj.Q + 
2 Yiy=iPrJ = 1 pi-j > (/ = 0, 1, . . . , c — 1), and is to be determined from the 
descendent distributions. It is noteworthy that the symmetry Pi{Si) = pi{—Si) on the 
boundary condition also restricts 7!'i{hi) to a symmetric function of the form of eq. (!33l) . 
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After the configurational average is performed, the cavity-field distribution Tii{hi) 
depends only on the distance, (?, from the boundary. Therefore, we hereafter denote 
7ii{hi) as Tigihi) and represent the distance of the origin from the boundary as g = L. 
The BP scheme assesses Pg+i using its descendents Pg. However, the only part relevant 
to the assessment of E{n) is that for pg-o, which is represented as 

P.+i;0 = ^ f ri, (34) 

l-2{l-e-^y) (i^j 

for the c = 3 case, being accompanied by an update of the partition function 

= |l - 2(1 - e-^y) (^^)'| , (35) 

and similarly for a general c. After evaluating pg.Q and Hg using this algorithm up to 
g = L — 1, the full partition function, S(?/), in the current limit n ^ and f] ^ oo 
keeping ?/ = n/? ~ 0(1) is finally assessed as 

m = ^L = ^UeJ^' 1 1 - 3(1 - e-'y) (^^f^) ' (1 + l>L_i;o) I . (36) 

For G C, eqs. f l34l) - fl36l) can be performed in a feasible computational time and 
therefore offer a useful scheme for examining RZs. This is the main result of the present 
paper. The concrete procedure to obtain RZs is summarized as follows: 

(i) To obtain a series of Pg.Q, eq. is recursively applied under the initial condition 
Po;o = until g reaches L—1. This can be symbolically performed by using computer 
algebra systems such as Mathematica. 

(ii) Using the series {Pgfl}, the moment Eg is recursively calculated by using eq. 
under the initial condition (l32l) until g becomes L — 1. Then, the full moment 
S(?/) = El is derived from eq. (|36ll using El-i and Pl~1;0- 

(iii) Solving El = with respect to a; = numerically. 

Although the right hand side of eq. flMl) is expressed as a rational function. El 
and E{y) are guaranteed to be certain polynomials of x since the contribution from the 
denominator is canceled in each step of eqs. (!35ll and (j36l) . The procedure of (i), (ii) 
and (iii) can be performed in a polynomial time with respect to the number of spins. 
However, for CTs the number of spins and the degree of the polynomial El increase 
exponentially as O [{{k — l)(c— 1))^) as L becomes larger, which makes it infeasible 
to solve El = for large L. For instance, it is computationally difficult to evaluate 
RZs beyond L = 7 for {k, c) = (2, 3) and L = 4 for {k, c) = (3, 4) by use of today's 
computers of reasonable performance. This prevents us from accurately examining the 
convergence of RZs to the real axis in the limit L ^ oo by means of numerical methods 
and analytical investigation for this purpose is non-trivial either. However, the data 
of small L still strongly indicate that the qualitative behavior of RZs can be classified 
distinctly depending on whether certain bifurcations, which are irrelevant to any RSB, 
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occur for the cavity field distribution in the hmit of L ^ oo. This imphes that RZs of 
the ladder and tree systems are related to no RSB. In the following sections, we give 
detailed discussions to lead this conclusion presenting plots of RZs. 



2.3. Remarks 

Before proceeding further, there are several issues to be noted. 



2.3.1. Uniqueness of the analytical continuation As already mentioned, analytical 
continuation from n G N to n G C cannot be determined uniquely in general systems. 
However, in the present system, we can show the uniqueness of the continuation. 
Therefore, the RS solution assumed above is correct. 

For this, let us consider the modified moment [(Ze~^^^)"]^/^, where A''^ is the total 
number of bonds. This quantity satisfies the inequality 



< 



< 



(Trl) 



Re(n) 



l/N 



2Re(n) ^ C'(e''l"l 



l/N 



(37) 



for finite A^. 
condition lib 



Suppose that we have an analytic function ipiji; N) that satisfies the 
[n; N)\ < ©(e'^'"'). Carlson's theorem guarantees that if the equality 



{n;N)-[{Ze- 



-l3NB\n1l/N\ 



holds for Vn G N, tp{n; N) is identical to [{Ze' 



-(3NB\n-ll/N 



for \fn G C. Because e~^^^ is a non- vanishing constant, this means that the 
analytic continuation of [Z"]^/^ is uniquely determined. This indicates that expressions 
analytically continued under the RS ansatz, namely eqs. (|30l) and (l36ll . are correct for 
finite A^ (or equivalently, finite L) although the analyticity may be broken on the real 
axis in the limit N oo. 



2.3.2. Relationship to other systems In addition to examining RZs for finite CTs and 
ladder systems, the relevance of RZs to the large system size limit will also be argued by 
comparison with known thermodynamic properties of relatives of CTs, namely Bethe 
lattices and regular random graphs. These are sometimes identified with CTs because 
the fixed point condition of the BP method is represented identically. However, we 
here strictly distinguish them. The definitions and properties of these systems are 
summarized as follows: 

• The Cayley tree (CT): A tree of finite size consisting of an origin and its neighbors. 
The first generation is built from c neighbors which are connected to the origin. 
Each site in the nth generation is connected to new c — 1 sites without overlap and 
all these new sites comprise the n + 1th generation. Iterating this procedure to the 
Lth generation, we obtain the CT, and the Lth generation becomes its boundary. 
For the boundary condition of eq. fl3T|) . which implies po;0 = in the expression 
of eq. ( IMl) . H(n) of this lattice is represented as a polynomial of x = e~^, which 
can be assessed by symbolic operations using eqs. (!34l) - (!36l) without evaluating the 
values of S(n). This property is very useful for investigating RZs. 
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• The Bethe lattice (BL): A lattice consisting of the first L' generations of a CT, for 
which L ^ oo is taken. Alternatively, we can define a BL as a finite CT of L' 
generation, the boundary condition of which is given by the convergent cavity field 
distribution of the infinite CT. Unlike for a CT, the boundary condition depends 
on y for a BL. Due to this difference, S(n) of this lattice cannot be represented as a 
polynomial and searching RZs becomes non-trivial. However, assessing the values 
of S(?T,) is still feasible computationally. 

• The regular random graph (RRG): A randomly generated graph under the 
constraint of a fixed connectivity c. Since there exist many cycles, assessing 

and RZs for this lattice is not feasible computationally for finite A^. In 
the limit N oo under appropriate conditions, however, it is considered that 
the RRG and the BL share many identical properties. Therefore, this lattice is 
sometimes identified with the BL and regarded as a solvable system [25l [26l [271 [28] . 
Nevertheless, we here distinguish between the two systems because the main 
purpose of this paper is to clarify the asymptotic properties of gN{n) from finite A^ 
to infinite A^, and our definition of the BL is useful to compare these limits. Here, 
the terminology "RRG" is used only to refer to systems of infinite size. 

2. 3. 3. Relevance of the boundary condition to the moment of the partition function The 
above mentioned distinction between the three relatives of CTs yields differences in the 
expression of the moment of the partition function, even while they share an identical 
cavity field distribution in the limit N —>■ oo. 

Equations ([MI) - ([36l) imply that gN^y) = N^'^ logH(y) for CTs is generally expressed 

as 

9^(y) = ^ - ^ - ^^9l'\y) + ^T.9.iy)^ (38) 

(ij) i M 

where g^^j'^iy) and gf\y) denote the contributions from the bond {ij) and the site i, 
respectively, and Cj is the number of bonds that site i has. The last term gfi{y) is the 
contribution due to the boundary fields. 

This is considered a generalization of a well-known property of free energies for 
cycle free graphs [29l [301 131] • For regular CTs, q = c holds if i is placed inside the tree, 
while Cj = 1 for the boundary sites. 

For a BL, the boundary condition given by the convergent solution of eq. flMj) . p^, 
which becomes a function of particularly simplifies the expression of eq. f[55]) as 

gf'{y) = rigi{y) + r^g^{y). (39) 

Here, n = (l + c(c - 2)-^ ((c - 1)^'"^ - l)) / (l + c(c - 2)-^ ((c - 1)^' - l)) and tb = 
1 — ri represent the fractions of the number of sites inside the tree and on the boundary, 
respectively, and 

gi{y) = lg'^'\y)-{c-i)g'^'\y), (40) 
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and 



c 



9B{y) = ^g^'\y) + g,iy), (41) 

represent contributions from a single site inside the tree and on the boundary. In general, 
g^'^\y) and g^^\y) are expressed as 

g^'\y) = log {Tr [p{S,r'piS2r'e''^'^'?'S] } , (42) 
g^'\y) = log {TrpiSy}, (43) 



where 



(2 cosh /?/;,) 

and Tr{h) is the distribution of the cavity bias, which is related to Tcih) as 



p{S) = I dh9{h) ^ (44) 



(45) 



TT{h) = / dhn{h) 



S (h- -^tanh~^(tanh/5Jtanh/3/i)^ 
For c = 3 in the limit f3n — > y, we have 

g^'\y) = loge^ (l - 1(1 - e-2^)(l - p,)^^ , (46) 

g^'\y) = log (^1 - ^(1 - e-'y){l-p^ni+p^)^ , (47) 

g,{y) =go- log (^1 - 1(1 - e-2^)(l - p^)'^ , (48) 

where go = log J dhP{h){2 cosh (3 h)^ is the contribution from a boundary spin and 
P{h) is the boundary-field distribution of the BL determined satisfying the condition 
7r(0) = P(0)/(/ P(/i)(2cosh/3/i)") =p,. 

Equation fl39|) represents a distinctive feature of cycle free graphs. In most 
systems, the contribution from the boundary becomes negligible as the system size 
tends to infinity. However, eq. (!39|) indicates that such a contribution does not 
vanish for a BL since re — > (c — 2)/(c — 1) remains of the order of unity even 
if = 1 + c(c — 2)~^ ((c — 1)^' — l) becomes infinite. Nevertheless, the complete 
separation of contributions between the inside and the boundary in this equation implies 
that it is physically plausible to use gi{y), instead of g^{y), in handling problems 
concerning the bulk part of the objective graph. Actually, such a replacement has been 
adopted in several studies on cycle free graphs [26l[32]. In general, gi{y) agree with g{y) 
of an RRG, which provides the basis of the correspondence between BLs and RRGs. 

In spin-glass problems on cycle free graphs, the replacement of g^{y) with gi{y) 
is crucial. To see this, we here investigate the large deviation properties of [Z^]. We 
denote the boundary condition as Psih) = IliGboundary '^ii.^i)- Equation ([6]) implies that 
S(?/) is expressed as S(|/) = / dhP-B{h) exp {—yEg{h)), where Eg{h) is the ground state 
energy when h is imposed on the boundary. For general systems, including a BL, this 
yields the identity 

y\d/dy) {y-'gN{y)) = N-'D{Pb\Pb) > 0, (49) 
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where PBih) = PBih) exp {-yEg{h)) /E{y) and D{Pb\Pb) is the Kullback-Leibler 
(KL) divergence between Pb(^) and PB{h). An imphcation of this relation from 
large deviation statistics is that the probability P{f) that Eg{h)/N equals / scales as 
P(/) ~ exp {NT,jq{f)) for large A^, where / and T^nU) are related by / = —{d/dy)gN{y) 
and SAf(/) = —y'^{d/dy) {y~^gN{y)) parameterized by y. The non-negativity of the KL 
divergence indicates that the rate function SAr(/) cannot be positive, which guarantees 
the normalization constraint j dfP{f) = 1. 

The constraint T^N^f) < is always satisfied even when oo. However, this is 
not necessarily the case when we take the thermodynamic limit limj^f^oo ONiv) = g{y) 
and then calculate the rate function as = —y'^{d/dy) {y~^g{y))- This function 

can be positive, and it can be shown that the condition = signals the 

onset of IRSB [121 ESj- The positive part of can be formally interpreted as the 
complexity or the configurational entropy of the metastable states for a single typical 
sample of couplings in the conventional IRSB framework [34j, as shown in fig. [2l In 
the IRSB framework, the critical condition = 0, which is alternatively expressed 

as {d/dy) {y^^g{y)) \y=ys = in general, corresponds to the typical state realized in 
equilibrium. 

The condition = has already been investigated for RRGs and indicates 

that IRSB transitions occur for some types of RRGs [261 [28]. However, it is considered 
that such a symmetry breaking cannot be detected by an investigation based on eqs. 
f lMl) - fl36l) because the boundary contribution is inevitably taken into account for a BL 
as well as for a CT. Actually, direct verification of < is possible for the c = 3 
case; details are shown in Appendix C This indicates that the possible RZs provided 
by the current scheme are irrelevant to IRSB. 




Figure 2. Schematic diagram of assessed using gi{y)- The rate function is 

continued to the complexity dX f — fs and the IRSB occurs at this point. The 
complexity vanishes at f = fd where the monotonicity of the free energy with respect 
to y = f3n breaks down. 
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Figure 3. Plot of RZs for a CT 

with c — 3. All the zeros lie on the 
line Ini(j/) — tt/2, as for a 2 x L 
ladder. 



Figure 4. RZs plot for a 3-CT 
with c — 3. A sequence of zeros 
approaches the real axis as the 
number of generations L increases. 
The arrow indicates the collision 
point expected from the study of 
the L oo limit in sec. 13.21 



3. Results 

3.1. Plots of the replica zeros for tree systems 

We here present the resuhs only for CTs. RZs plots for ladder systems are summarized 



m 



Appendix A The plots for a CT and for a 3-CT with c = 3 are shown in figs. [3] and 
m respectively. Note that IRSB occurs in RRGs with the same parameters. The critical 
values are ys = 0.41741 and oo for the RRG counterparts of a CT and 3-CT with c = 3, 
respectively. 

Figure [3] shows that RZs of the c = 3 CT lie on a line Im(?/) = tt/2. Interestingly, 
this behavior is the same as the 2 x L ladder case, the plot of which is given in 



Appendix A This result indicates that there is no phase transition or breaking of 
analyticity of g{n) with respect to real y. This is in accordance with the argument on 
the boundary contribution mentioned in the previous section. 

On the other hand, for the 3-CT case in fig. |U a sequence of RZs approaches a point 
yc on the real axis from the line Im(y) = tt/2 as the number of generations L increases, 
although the value of yc is far from y,^ = oo. A similar tendency is also observed for a 
CT and 3-CT with c = 4, plots of which are presented in figs. [5] and [6l respectively. 
The IRSB critical values are = 0.38926 for the CT and ys = 1.41152 for the 3-CT. 
Again, these values are far from the values of yc, which can be observed in figs. [5] and [61 

These results indicate that certain phase transitions occur for some CTs, although 
they are irrelevant to IRSB. It is difficult to identify the critical value yc from the plots 
because of the computational limits. Instead, in the following subsection we investigate 
the L —>■ oo limit of these models. The arrows in figs. SHE] represent the transition 
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Figure 5. RZs of a CT with 
c = 4. We consider only an 
L-generation branch in this case 
because of computational limits. 
RZs approach the real axis as L 
increases around tJc ~ 0.5. The 
arrow indicates the location of 
the singularity of the cavity-field 
distribution. 



Figure 6. RZs of a 3-CT with 
c — 4. We consider only an 
L-generation branch. The zeros 
approach the real axis around yc ~ 
1.1. There are two singular points 
of the cavity field distribution 
in this case, both of which are 
indicated by arrows. 



points He determined by this investigation. 
3.2. Phase transition on the boundary of a BL 

In order to identify the value of Uc, we take the hmit L — oo by equating Pg+i-o and Pg.^Q 
in the iterative equation of Pg-o, which yields the boundary condition p^ of the BL. For 
a c = 3 3-CT, the iterative equation is given by 

_ [Plfi + 2P,;0(1 - Pg;0)Y + je-^^l " P.flY 

" l-|(l-P.o)ni-e-^) • 

A return map of the recursion of Pg-o and the convergent solution p^ are presented in 
figs. [7] and [HI respectively. The return map shows that there are three fixed points for 
X ^ 2.35, while p = 1 is the only fixed point for x ^ 2.35. This situation is in contrast 
to the c = 3 CT case, in which the cavity- field distribution uniformly converges to an 
analytic function: 



ix'^ 



P* = — ;:2^ ' (51) 



2 + - /FT 
2{l-x\ 

which can be derived from eq. This implies that when eq. (!3T!) is put on the 

boundary of the CT, the boundary condition of the BL, which was obtained by an 
infinite number of recursions L — L' oo, exhibits a discontinuous transition from 
< 1 to = 1 at a; ~ 2.35 ?/c ~ 0.85 as y is reduced from the above. Actually, in 
fig. m RZs of the c = 3 3-CT seem to approach yc ~ 0.85, marked by an arrow. This 
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0.2 0.4 0.6 



Figure 7. Return map of a 3-CT 
with c = 3. Tlie convergent point 
of tlie recursion discontinuously 
changes depending on x. The sohd 
hne represents the function f(jp) = 
P- 




\ \ \ 



\ 



\\\ 
\ \ \ 
\ \ \ 

V\ 



Figure 8. Asymptotic behavior 
of pgfi of a 3-CT with c — i. A 
finite jump of p occurs at a; w 2.35. 
The solid hne denotes the h ^ oo 
solution . 




Figure 9. Return map of a CT 
with c = 4. The stable fixed point 
is unique but shows a singularity at 



Figure 10. Asymptotic behavior 
of Pg-fi of a CT with c = 4. 
In the thermodynamic limit, Pgfi 
is continuous but the derivative 
becomes discontinuous at x — \/3- 



indicates that RZs obtained by our framework are relevant to the phase transition of 
the boundary of a BL, which is not related to IRSB. 

The same analysis for a c = 4 CT shows that bifurcation of another type can occur 
for even c. For this model, the recursive equation of Pg-o has a trivial solution = 
for Vx, which is always the case when c — 1 is odd. The return map and plots of are 
shown in figs. M and [TOl respectively. 

These figures indicate that there exists a continuous transition from = to 
p^, > at a certain value of x, which can be assessed as Xc = =^ yc ^ 0.5. This is 
consistent with a certain sequence of RZs approaching the real axis around yc ~ 0.5 in 
fig. which supports the analytical assessment of the critical points. 
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In general, the discontinuous transition appears for cases of A; > 3 spin interactions 
and the continuous transition occurs when c is even. Actually, for a c = 4 3-CT, 
both discontinuous and continuous transitions occur at x ~ 0.86 =^ y ^ —0.15 and 
X = 3 ^ y ^ 1.1, respectively. Figure [6] shows a sequence of RZs approaching y ^ 1.1, 
while it is difficult to clearly identify a sequence converging to the other critical point 
y ^ —0.15. We consider that this is because the system size is not large enough, since 
a portion of the RZs in the left shows a tendency to approach the real axis, though 
further increase of the system size is practically unfeasible due to the limitations of 
current computational resources. 

In conclusion, the analysis shown in this section indicates that RZs of CTs are 
related to the phase transitions on the boundary of a BL. Regardless of the type of 
transition, a sequence of RZs approaches a critical point on the real axis when the BL 
provided from a CT in the limit L ^ oo exhibits a phase transition on the boundary. 



4. Discussion 



4.I. Possible link to AT instability 



The AT condition, which is critical for FRSB, has not yet been characterized for sparsely 
connected spin models. In fact, previous research has found that critical values of the 
continuous transitions from = to > are candidates for those of the AT condition 
for systems of even c [28]. This motivates us to further explore a possible link between 
RZs and the AT instability. 

Divergence of the spin-glass susceptibility of the root site is often adopted as 
the critical condition of the AT instability for BLs [351 ISSl [371 IM]- Generalizing the 
condition to the case of finite n, we obtain 



XSG 



dh, 



(52) 



where [(■ ■ ■)]^ means an average with respect to a modified distribution of coupling and 
boundary field 

F({J.,},{/..})Z-({J,,},{/i,}) 



Y.{J^,}P{{J^^},{h^})Z-{{.h,},{K}y 



(53) 



This definition is reasonable because eq. (]52l) correctly reproduces the AT condition of 
fully connected systems for finite n in the limit of infinite connectivity c — > cxo [33| [39] . 

In a cycle-free graph, an arbitrary pair of nodes is connected by a single path. Let 
us assign node indices from the origin of the graph to a node of distance G along the 
path as (7 = 1, 2, . . . , G. For a fixed set of couplings and boundary fields, the chain rule 
of the derivative operation indicates that 



d {So) d {So) dho dho dh 



d {So) dho yr dhg 
dho dh A A 



G 



dh, 



G 



dhn 



dho dh. 



dh 



G 



I dh„ 



g=l 



dh„ dh„ 
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d (So) -Q dhg^i 



dh. 



5=1 



(54) 



as hg depends linearly on hg as hg = hg + rg, where represents a sum of the cavity 
biases from other branches that flow into node g. For a BL of {k,c) = (2,3), the BP 
update yields an evolution equation of the cavity bias 

1 



h 



9-1 



/3 



tanh ( tanh(/3Jg) tanh{/3{hg + rg)) 



oo 



(55) 



sgn i^Jg{hg + rg)j ( \hg + rg\ > 1 

Jg{hg + Tg) ( otherwise ) 

where Jg denotes the coupling between nodes g — 1 and g, and similarly for other cases. 

To assess eq. ( l52l) . we take an average of the square of eq. ( I55l) with respect to 
the modified distribution Pn{{Jij},{hi}). Here, Vg can be regarded as a sample of a 



stationary distribution determined by the convergent solution of eq. (IMI) for the BL. As 



is limited to being an integer and | J, 




1, eq. fl55l) gives 

( fig + Tgl > 1 ) 

{\h9 + rg\ = l) , 
( otherwise ) 



(56) 



where the value of or 1 for the case of \hg + rg 
value of hg. When hr, 



1 is determined depending on the 
(and \rg\ = 1), the values and 1 are chosen with equal 
probability 1/2 since the sign of the infinitesimal fluctuation of hg, 6hg, is determined 
in an unbiased manner due to the mirror symmetry of the distribution of couplings. On 



the other hand, under the condition of n^=q+i dh^^i/dhk ^ 0, the case of \h 



9+1 



(and rg = 0) always yields 
for \hn\ 



dhg_i/dh. 



■ 1. This is because hg5hg < is guaranteed 

1 under this condition. 
Equation (156!) indicates that the assessment of eq. (15^ is analogous to an analysis 
of a random-walk which is bounded by absorbing walls. We denote by P(g^o) the 



probability that 



dhg^i/dh^ 



never vanishes during the walk from G to and the value 



of n^=i \9hg_i/dhg\ is kept to unity. This indicates that 



d{So) 
dh. 



G 



OC P, 



(57) 



holds. Summing all contributions up to the boundary of the BL yields the expression 

L' 



X5GCX^(A:-l)«(c-l)«P( 



G=0 



The critical condition for convergence of eq. fl58|) in the limit L' — > oo is 
log {{k - l)(c - 1)) + hm i logP(G_o) = 0. 

G^OO (jr 

This serves as the "AT" condition in the current framework. 



(58) 



(59) 
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For a BL, eq. (1591) can be assessed by analyzing the random walk problem of eq. 



( !56|l . as shown in Appendix D We evaluated the critical ijat values of eq. fl59l) for several 
pairs of {k,c), shown in Table [H along with other critical values. These results show 



(fc,c) 


Vat 


Vc 


Vs 


(2,3) 


0.54397 


none 


0.41741 


(2,4) 


0.89588 


log 0.54931 


0.38926 


(3,3) 


1.51641 


0.85545 


oo 


(3,4) 


1.35403 


-0.15082, log3 ^ 1.09861 


1.41152 



Table 1. Relevant values of y. Note that each kind of y is calculated using different 
models. The IRSB transition point is for RRGs and yAT is for RRGs or BLs. The 
singularity of the cavity-field distribution y^ is common for all the models. 



that the values of t/c, which signal the phase transitions of the boundary condition of the 
BL, agree with neither yj^j' or i/s, implying irrelevance of RZs to the replica symmetry 
breaking. 

The irrelevance of RZs to the AT instability may be interpreted as follows. We 
can link the spin-glass susceptibility to gN{n) in general by considering the following 
extension: 

m N \\ 

Tr exp I -/? ^ H{Sn + E E ^i^i i^' e-^^^-"^ , (60) 

a=l 1=1 a<b / / 

by breaking the replica symmetry introducing replica symmetric interactions among 
m out of n replica systems with coupling F = (Fi,F2, . . . ,-Fjv). Obviously, gN^n) = 
'gN{,F = 0;m,n) and gniji) = gN{F;m = l,n) hold. Analytically continuing eq. (!60ll 
to n, m G M and expanding the obtained expression around F = for m ~ 1 yields 

m — 1 




NgNiF; m, n) ^ NgNin) + 



-F XsgF + (higher orders). 



(61) 



where xsg = ([((5';S'fc) — {Si) (^fc))^]^) represents the spin-glass susceptibility matrix. 

Equation fl6T|) implies that the divergence of the spin-glass susceptibility is linked 
to analytical singularities of \im.M -,00!] N{F\m,n) for m 7^ 1. However, for m = 1, 
which corresponds to gN^n) examined in the present paper, it is difficult to detect the 
singularity because the factor m — 1 with F'^xsgF makes the divergence of the spin- 
glass susceptibility irrelevant to the analyticity breaking of g{n) = limAr^oo 5'7v(^)- A 
possible solution is to consider systems of m 7^ 1 in the framework of IRSB. However, 
an examination along this direction is beyond the scope of the present paper. 



4-2. Physical implications of the obtained solutions 

We concluded that bifurcations of the fixed point solutions of the BP update correspond 
to phase transitions of the boundary condition of a BL and are not relevant to either 
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IRSB or FRSB. Before closing this section, we discuss the physical implications of the 
obtained solutions. 

A naive consideration finds that the solution of Poo;0 = p = 1 corresponds to 
a paramagnetic phase implying that any cavity fields vanish and therefore all spin 
configurations are equally generated. Note that this phase is of the ground states in the 
limit /3 — > oo and is different from the usual temperature-induced phase. 

For finite p < 1, relevant fractions of the spins can take any direction without 
energy cost because the cavity field on the site is 0. This imphes that the ground state 
energy is highly degenerate, which means that this solution describes a RS spin-glass 
phase. Actually, it is easy to confirm that the following equahty holds: 

= = Tr S^S^^p.iS) = 1 - p,,o. (62) 

Hence, the singularity of the cavity- field distribution in the limit — > oo can be regarded 

as the transition of the spin- glass order-parameter. A finite jump of TToo{h) for the k = 3 
case is the first-order transition from the RS spin glass to paramagnetic phases, and 
such a transition is also observed in the mean-field models. The transitions from p = 
to finite-p for the c = 4 case can be regarded as a saturation of q to q^A = 1- We infer 
that these are the transitions from RS to RS phases. Notice that such a transition has 
not been observed for infinite-range models. Our results indicate that this q — I phase 
appears only when c is even. This means that such a phase is highly sensitive to the 
geometry of the objective lattice. This may be a reason why such a transition has not 
been observed in other models. 

5. Summary 

In summary, we have investigated RZs for CTs and ladders in the limit T, n — 0, /9n — > 
y ~ 0(1). Most of the zeros exist near the line Im(y) = 7r/2 in all cases investigated; in 
particular, for the (/c, c) = (2, 3) CT and the width-2 ladder all the zeros he on this line. 
For the width-2 ladder we have proved that the free energy is analytic with respect to y 
in this model. On the other hand, for some CTs, a relevant fraction of the RZs spreads 
away from the line lm{y) = tt /2 and approaches the real axis as the generation number L 
grows. This implies that g{n) has a singularity at a finite real y in the thermodynamic 
limit. A naive observation finds that the RZs collision points correspond to phase 
transitions of the boundary condition of the BL. We have compared them with known 
critical conditions of IRSB and FRSB and concluded that these conditions are irrelevant 
to the behavior of RZs. This is consistent with the absence of RSB in CTs reported in 
some earlier studies. 

To fully understand and use the replica method, as well as mathematical 
verification, an description of the physical significance of the method is required. We 
hope that our results presented in this paper lead to a deeper understanding of the 
mysteries of the replica method. 



Complex Replica Zeros of ± J Ising Spin Glass at Zero Temperature 



21 




i-1 



Figure Al. Unit cell of a 2 x L ladder. 
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Appendix A. Results for ladder systems 

We first explain the procedure to obtain the RZs of ladder systems. For a 2 x L ladder 
system, the BP equation can be derived in a similar manner to the CT case. We trace 
out the two spins of the previous generation to the iih generation, as in fig. lAll This 
yields an expression corresponding to eq. (IHl) as 



q1,2 



= AeM(3{r + Ji.i)SlSj). (A.l) 

From simple algebra we obtain 
^ 1 

J._i = -tanh"^(tanh/3irnanh/3i^'^tanh/3Ji_i), (A.2) 
/-^ 

cosh.B J i^i cosh. BK^ cosh. BK"^ , 

A = 4 ^ . (A.3) 

cosh/?Jj_i 

This shows that the effective bond Jj between 5*/ and Sf becomes 

J^ = r + X^l- (A.4) 

These relations indicate that the one-site marginal distribution pi for trees is replaced by 
the two-site marginal distribution Pi{S}, Sf) for a width-2 ladder. From the symmetries 
of the original model, we can specify the form of this distribution as 
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This expression can be interpreted as showing that the effective bond ffuctuates by 
quenched randomness. In a similar way to the tree case, the iterative equation for vr(J) 
is derived as 



ni{Ji) oc / dJi_iHi_i{J, 



5{Ji — J* — Jt 



i-i 



cosh pji 
cosh /9 J, 



(A.6) 



and that for the effective partition function is 



dJi^i'Ki^i{Ji^i 



{2 cosh f3KT {2 cosh /3K 



2\n 



cosh 13 Ji 



(A.7) 



cosh /?Ji_i 

In the limit (3 —>■ oo, (3n y, we can derive the following formulas from above equations: 

1 - Pg-fi 



-9+1 



o3y 



.0 + 2(1 



)(l-P9;o) 



(A.8) 
(A.9) 



Using these relations, we can symbolically calculate the as a polynomial of x = 
and obtain the RZs by numerically solving = as the CT cases. 

On the other hand, for larger-width ladders, the number of spins added by an 
iteration is greater than 2 and many-body interactions appear. It makes the problem 
complicated and simple relations like flA.4l) cannot be obtained. Hence, when treating 
larger-width ladders, we directly use the BP formula for a given sample {Jij} to obtain 
the ground state energy [IQ] and numerically assess the distribution of the ground state 
energies P{Eg) by enumerating all the configurations. Using the distribution P{Eg), the 
RZs equation is derived as 

^ P(Eg)e-^^« = 0. (A. 10) 

Eg 

This equation is solved numerically in the same way as the other cases. Note that 
computational times required in the counting process to obtain P{Eg) exponentially 
increases as L grows, which makes it infeasible to obtain for large L. 

Next, we present the plots of RZs for ladder systems with brief discussions. Figure 
IA2I shows the plot for a 2 x L ladder with the boundary condition pq.q = and Eq = x. 
Notice that all RZs lie on a line Im(?/) = tt/2. This fact can be mathematically proven, as 



detailed in [Appendix B[ The physical significance of this behavior is that the generating 
function gN{n) is analytic with respect to real y even for the N —>■ 00 limit. We have also 
investigated a 3 x L ladder and found qualitatively similar results as for the width-2 case. 
For a width-4 ladder, the RZ plot is given in fig. IA3[ We can observe that some zeros 
approach the real axis around Re(?/) ~ 1.2, but the rate of approach decreases rapidly 
as L grows. This implies that the RZs do not reach the real axis, which agrees with a 
naive speculation that ladders are essentially one-dimensional systems and therefore do 
not involve any phase transitions as long as the width is kept finite. 
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Figure A2. RZs for ladders 
with 2 X L. All the zeros lie on 
Ini(y) = n/2 and never reach the 
real axis of y = (3n. The inequality 
R-e(?/) < log 2a/2 holds, as shown in 
[Appendix B[ 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Re[y] 

Figure A3. Zeros of width- 
4 ladders. Some of the zeros 
approach the real axis around 
Re(t/) w 1.2, but the rate of 
approach rapidly decreases as L 
grows. 



Appendix B. Location of replica zeros of a width-2 ladder 

We prove that all RZs of a 2 x L ladder lie on the line Im(?/) = 7r/2 for any L. We 
introduce the notation 

/ N ni(x) ,^ ^ , 

di{x) 

where di and ni are polynomials of x = and ni{x)/di{x) is assumed to be irreducible. 
The outline of the proof is as follows. First we present the general solution of pi and show 
that the denominator di has 2F{{1 + l)/2) roots which are all purely imaginary. The 
function F{1) is the floor function, which is defined to return the maximum integer i in 
the range i <l. Also, we show that the number of nontrivial solutions of = is equal 
to 2F((/ + l)/2) and can be factorized as Ci{x)di{x), where Ci{x) is a polynomial of 
X. From the correspondence of the numbers of the roots, we conclude that all the zeros 
of 2/ are equivalent to the roots of di{x) and Ci{x) takes the form ax^. 

The iteration flA.81) for pi has a solvable form and its general solution is given by 



p, = ^-M^f) ^ , (B.2) 



4'(2 + x2 - + 8) - (2 + x2 + xVa;2 + 8)/i(x, 



where 



h{x) = -A- x{x + V^^) = 4 ^ ^ v^g^ _ ^g_3^ 

X — V + 8 

The roots of the numerator in eq. (16.20 can be easily calculated as 

±2v/2, ('72'" + l), (B.4) 

^ 0, ±2^2^ (/ = 2m) ' ^ ' 
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where i denotes the imaginary unit and m is a natural number. Then, we concentrate on 
finding the roots of the denominator in eq. flB.2p except for those of the numerator flB.4p . 
From numerical observations in sec. 3, we found that any of the roots x*, which satisfy 
= 0, are purely imaginary and bounded by < 2a/2. Hence, we assume these 
conditions and perform the variable transformation z = —xi. Equating the denominator 
of eq. (IB. 21) to 0, we get 



hi-ixW fVS^ + izV l-z^-i^f^^ 



4 ) \^/s^^-izJ 2-z^ + i^/. 



y2 



z 



2 



(B.5) 



We now enumerate the number of solutions under conditions that z is real and bounded 
as — 2a/2 < z < 2\/2. Under these conditions, we can transform eq. (IB. 50 into a simple 
form by using the polar representation. The result is 



where VS - z'^+iz = rie*^i (-vr < 6'i < vr) andr2e*^2 = 2-z^-iV8 - z"^ (-tt < 02 < tt). 
While z varies from —2^/2 to 2\/2 continuously, the radius ri stays at a constant 2^2 
and the argument 6i varies from —tt/2 to 7r/2 in the positive direction. In the same 
situation, 62 changes from +7r to —it in the negative direction. The radius r2 is not 
constant, but is finite in this range. The variables 61 and 62 are obviously continuous 
and monotonic functions of z. Therefore, the argument of the left-hand side of eq. (IB.6P 
starts from 6 = and rotates with angle 2/7r in the positive direction and the counterpart 
of the right-hand side varies from the same point 6 = to —Att. This means that there 
are / + 1 values of z where the factor {26i{z) — tt)/ becomes equal to 262{z) except for 
trivial solutions z = ±2\/2. When / is even, these solutions contain a trivial solution 
z = 0, which can also be confirmed from eq. (IB. 21) . Hence, the number of nontrivial 
roots of di becomes / + 1 for odd / and / for even /, which is equivalent to 2F((/ + l)/2). 

As already noted, the number of nontrivial solutions of = is equal to 
2F{{1 + l)/2). This can be understood by considering that the number of terms of 
[Z^] is determined by the maximum number of defects n^. In the 2 x / ladder case, the 
value of is given by F{{1 + l)/2) and the number of terms is + 1. The highest 
degree of the relevant polynomials for RZs comes from the difference between the highest 
and lowest ground-state energies and is given by 2nd = 2-F((Z + l)/2), which yields the 
number of nontrivial solutions of S/ = 0. 

Finally, we prove that takes the form Aix'''di{x) by induction. From eqs. (lA.SP 
and (IA.9|) with the initial conditions po;0 = 0, So = x, we derive 

Pi = ^^' 2i = ^x2(x2 + l), (B.7) 

which satisfies the desired form. Assuming that the condition S/ = Ax'''di{x) is true for 
I = k, v/e substitute this expression into eq. (1A.9P to get 

= ^Ax'^+' { {x' - l)nk + (1 + x^)dk} . (B.8) 
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Equation (lA.Sp can be written as 

P'+^ - (x2-iK + (i + x2)4 " 477' ^ 

which gives 

(x^ - l)nk + (1 + a;^)4 = Cfc+i(x)4+i(a;), (B.IO) 

where Ck+i is a polynomial and satisfies c^+i = {dk—rik) / ri^+i- Substituting this relation, 
we can rewrite eq. fIB.Sp as 

Hfc+i = ]^Ax^*^+\k+i{x)dk+i{x). (B.ll) 

As we have already shown, the number of nontrivial zeros of H^+i is equal to that of 
4+1- This means that Ck+i cannot have nontrivial roots and hence Ck+i takes the form 
Ax^. This completes the proof by induction and demonstrates our proposition that all 
RZs for a 2 X L ladder have a constant imaginary part i7r/2. 

Appendix C. Rate function for a CT with c = 3 

We here calculate the generating function giiv) for finite L. Consider an L-generation 
branch of a c = 3 CT. An explicit form giill) is easily derived from eq. fl35|) as 

^^(a;) = ^"^^^ (l-2-^)logx+ L-_y !^,(C.l) 

i'^v / 2 — \ 2 + — 1 / o 4 — 2^ + / J 2* 

i=0 

where a; = and 

5(0 = log Ho, /i = /i(x,pi;o) = 1 - -(1 - -]9^;o)^ (C.2) 

using the same notations as in sec. [H The rate function with finite generations L is 
given by 

dgo 



S.(a:) = ^,;^U-a:logx.^^ 



L-l 



4 _ 2-i+i Z-j 



2 = 

where the factor Ci{x) is given by 

Q{x) = ^^ + ^ = {1- x-'){l - p,o)^ - ^=^(1 - P.o)^. (C.4) 
apj.Q ax ox ax 

Let us denote Soo(a^) = liiiii_^oo ^^(a;). Because the inequality Sqo < always 
holds, the IRSB transition does not occur as long as the condition T,oc{x) = is 
satisfied. 

In the range y > 1 < x, the factor /j is bounded as 1/2 < /j < 1. This 
guarantees the uniform convergence of gL^x). The boundedness of (dpi-o/dx) can also 
be shown with some calculations. These conditions guarantee that t converges to 

a function Eqo uniformly. Hence, from elementary calculus, the equality = Soo(x) 
holds, which implies the absence of IRSB. The same conclusion is more explicitly derived 
for a BL because fi does not depend on i. 
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Appendix D. AT condition for the {k,c) = (2, 3) case 

We here evaluate the AT condition for a BL with {k,c) = (2,3). To evaluate P(g->o), 
we construct the transition matrix of our random- walk problem. For a given [hg, hg^i), 
the posterior distribution of Vg is given as 

p{rg\hg) = p{rg,hg)/p(hg) oc d ^ I ^''^ I ^ 1^ ( r, ) , (D.l) 

where p{rg) is the prior distribution of Vg. This enables us to derive the concrete 
expression of p{rg\hg), summarized in Table IDll We can distinguish three states of 



rg \ hg 


1 





-1 


1 




1 -P6 


{i-Pb)e-'y 


(1 + Pb) + (1 - Pb)e-^y 


2 


(1 + Pb) + (l - Pb)e-^y 





'2pt 


Pb 


'2pb 


{1 + Pb) + (i - Pb)e-^y 


{I + Pb) + (i - Pb)e-^y 


-1 


{i-Pb)e--'y 


I-Pb 


1 -Pfe 


{I + Pb) + {I - Pb)e~^y 


2 


{1 + Pb) + (i - Pb)e-^y 



Table Dl. Values of p{rg\hg) for (fc, c) = (2, 3). The symbol pi, is the probability that 
the cavity bias takes the value 0. 

the walker at the g-step as follows: 
|1): The walker has already vanished. 
|2): 
|3): 

Hence, using the relation fl56|) . the transition matrix T can be written as 



The walker survives and \hg\ = 1. 
The walker survives and \hg\ = 0. 



/ 1 Pi,i |Pi,o X 2 \ 



(D.2) 



Po,i ^Pi,o X 2 

where p r represents p(ro\ho) and the condition p r = p t applies. When 

^ ^Si"9 i ^ J \ yi M/ J rg,ng ■> —rg,—ng 

\hg + rg\ = 1 and \hg\ = 0, the states |1) and |2) occur with equal probability 1/2, 
while 1 2) is always chosen when \hg + rg\ = 1 and \hg\ = 1 as exlained in sec. 14.11 
This matrix has three eigenvalues: Ai = 1, A2, and A3. The eigenvector of the largest 
eigenvalue Ai = 1 corresponds to the state |1) or the vanishing state. Hence, the 
surviving probability P(g-^o) is given by 1 — (1|G), where \G) is the state of the walker 
at the G step. For large G, the relevant state is of the second-largest eigenvalue A2, and 
we get 

PiG^o) ~ A^. (D.3) 

Using the stationary solution flSTj) . we obtain P(g^o) as a function of x = e^. The AT 
condition becomes 

XSG OC - l)''(c - l)''P(G^o) ^ 00 ^ (fc - l)(c - 1)A2 > 1. (D.4) 

G 
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This condition is easily examined numerically and we can verify that the AT instability 
occurs at yAT ~ 0.54397 for (fc, c) = (2, 3). 
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